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By introducing Dirac (5-function in superhigh magnetic field, we deduce a general formula for pressure of degenerate and 
relativistic electrons, Pe, which is suitable for superhigh magnetic fields, discuss the quantization of Landau levels of 
electrons, and consider the quantum electrodynamic(QED) effects on the equations of states (EOSs) for different matter 
systems. The main conclusions are as follows: the stronger the magnetic field strength, the higher the electron pressure 
becomes; compared with a common radio pulsar, a magnetar could be a more compact oblate spheroid-like deformed 
neutron star due to the anisotropic total pressure; and an increase in the maximum mass of a magnetar is expected because 
of the positive contribution of the magnetic field energy to the EOS of the star. Since this is an original work in which 
some uncertainties could exist, to further modify and perfect our theory model should be considered in our future studies. 
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1 Introduction 

Pulsars are among the most mysterious objects in the uni¬ 
verse that provide natural laboratory for investigating the 
nature of matter under extreme conditions, and are univer¬ 
sally recognized as normal neutron stars (NSs), but some¬ 
times have been argued to be quark stars (Du et al. 2009, 
Lai et al. 2013, Xu et al. 2013). The equation of state (EoS) 
of matter under exotic conditions is an important tool for 
understanding of the nuclear force and for astrophysical ap¬ 
plications. The Fermi energy of relativistic electrons Ep{e) 
is one of most important and indispensable physical param¬ 
eters in EoS, and affects direct weak-interaction processes 
including modified URCA reactions, electron capture (e.g., 
Gao et al. 2011a, 2011b, 2011c, 201 Id, 2012a, 2012b; Liu 
2012, 2013, 2014, 2015; Du et al. 2014). These influences 
will change intrinsic EoS, interior structure and heat evolu¬ 
tion, and even affect the whole properties of the star. 

As we know, for degenerate and relativistic electrons 
in /3—equilibrium, the distribution function f{Ee) obeys 
Fermi-Dirac statistics: f{Ee) = \/[Exp{{Ee —^ie)lkT)-\- 
1), k represents Boltzmann’s constant, and is the electron 
chemical potential. If T —0, pe is also called the electron 
Fermi energy, Ep{e), which presents the energy of high- 
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est occupied states for electrons. The electron Fermi energy 
Ep{e) has the simple form 

Ep{e) = {pl{e)c'^ + , (1) 

with pp{e) being the electron Fermi momentum. 

In the context of general relativity principle, the matter 
density is defined as; p = ejc^, the total energy density, 
including the rest-mass energies of particles. Using the basic 
thermodynamics, we obtain the relation of the total matter 
pressure P and matter density p in a common NS, 

2 d{e/nB) 

P{nB) = nB—: - , 

driB 

p{nB) = e{nB)lc^, => P = P{p) . (2) 

From the above equation, it is obvious that P solely depends 
on p. Theoretically, we can obtain the value of Ep{e) by 
solving EOS in a specific matter model. The pressure of de¬ 
generate and relativistic electrons, Pe, is another important 
and indispensable physical parameter in EoSs of a NS. Pe is 
one of important dynamical pressures against a NS’s gravi¬ 
tational collapse, and affects the structures and properties of 
the star, substantially. 

Thompson and Duncan (1996) predicted that superhigh 
magnetic fields (MFs) could exist in the interiors of magne- 
tars with a typical surface dipolar MF, B ^ 10^^ to 10^® G 
(Thompson & Duncan 1996). Superhigh MFs have effects 
on EoSs of a NS, as well as on its spin-down evolution (e.g., 
Gao et al. 2014, 2015). Recently, Franzon et al. (2015) stu- 
did the effects of strong MFs on hybrid stars by using a full 
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general-relativity approach, and pointed that the MF could 
cause the stellar central density to be reduced, inducing ma¬ 
jor changes in the populated degrees of freedom and, po¬ 
tentially, converting a hybrid star into a hadronic star. In ac¬ 
cordance with the popular point of view, the stronger the 
MF strength, the lower the electron pressure becomes. With 
respect to this view, we cannot directly verify it by exper¬ 
iment in actual existence, owing to lack of such high MFs 
on the earth. After a careful check, we found that popular 
methods of calculating Ep (e) the electron Fermi energy are 
contradictory to the quantization of electron Landau levels. 
In an extremely strong MF, the Landau column becomes a 
very long narrow cylinder along MF, If we consider Dirac 
(5-function in superhigh MFs, all the results should be re¬ 
considered. 

In Sec. 2, we deduce an equation of Pg in superhigh 
MFs; in Sec. 3, we consider QED effects on EOSs of NS 
matter, and discuss an anisotropy of the total pressure; In 
Sec. 4 we discuss our future work on improving our model, 
and present conclusions in Sec. 5. 


2 Deduction of the pressure of electrons in 
superhigh MFs 


The relativistic Dirac-Equation for the electrons in a uni¬ 
form external magnetic field along the z—axis gives the 
electron energy level 


Ee 


\mtc^ 


2P 

(1 + 

-Dcr 


) 



(3) 


where i^ = n+^+a is the quantum number, n the Landau 
level number, cr = ± ^ the spin quantum number (Canuto 
& Ventura 1977), and pz is the z-component of electron 
momentum, and may be treated as a continuous function. 
Combining Per = mlc^/eh with = ehl^vtieC gives 

Pg = TOgC"* +FzC^ + ‘2i'2meC^ PgB 

_ 2 4 I 2 2 I 22 /A\ 

= m^c +p^c +pj_c , (4) 


where is the magnetic moment of an electron, and 
p± = mgc(2r/P*) 2 ). The maximum electron Landau level 
number Umax is uniquely determined by the condition 
[P'f{z)c]‘^ > 0 (Lai & Shapiro 1991, Gao et al. 2013), where 
ppiz) is the Eermi momentum along the z—axis. The ex¬ 
pression for iZmax can be expressed as 

^max{^ — 2 ) 




1 r.PF(e). 


2B* ruec^ 


UleC 2 2 


1 


^Epie 


= - 1 - (^f] 

2B* rrieC^ rrieC 

^max — 2 ) 

^2B*'-^ ruec^ 

- T 

^ 2P* ^ TOgC^ 


- 1 - 


r]-i + ^+ J] 

TTleC 2 2 


l-( —f] 

TTleC 


(5) 


According to the definition of pF(e) in Eq.(l), we ob¬ 
tain pF(e) = pF{e)c if electrons are super-relativistic 
{Ep{e) me(?). In the presence of a superhigh ME, 
B ^ Per, pF(e) nieC^), we have 

^maxi^ 2^ 2^ 


~ 7nf[- 


1 




(6) 

The maximum of p± for electrons in a superhigh ME is 
pl{Tnax)c^ = 2u'^^^2meC^PeB , (7) 

where the relation of 2p^Bcr/meC^ = 1 is used. Inserting 
Eq.(6) into Eq.(7) gives 

1 .PF(e) 

2B* ^ mgc2 


p‘]_{max)c^ = 2 X x 2meC^PeB 


~ Per X X 2mgC^ 


mlc^ 

eh 


eh 

UleC^ ' 2meC 

,Ep(e).^ „ e, eh 

< X 2mec^- - 

rrieC^ IrUeC 


= El{e) .(8) 


In superhigh MEs, Ep{e) is determined by 

PF(e):x 43.44(A)i/4(^_^)l MeV , (9) 

Bcr Po U.UOOO 

where pQ = 2.8xl0^^g cm^ is the standard nuclear density 

(Gao et al. 2012b). Thus, we obtain 

/ \ / \ ^f(6) 

p±[max) = pp(e) ~- 

c 

= 43.44 X A)iMeV/c (P* > 1) (10) 


'0.0535 Po Be 

As pointed out above, when n = 0, the electron Landau 
level is non-degenerate, andp^ has its maximum pz(TOaa;), 

i ^ ^ ^ pF(e) 

Pzimax) =pp(e) ~ - 

c 

= 43.44 X A)lMeV/c (P* > 1) (11) 


'0.0535 Po Be 

Erom Eq.(9) and Eq.(lO), it’s obvious that pz{max) = 
p±{max) = pp{e). The reason for this is that in the in¬ 
terior of a magnetar, electrons are degenerate and super- 
relativistic, and can be approximately treated as an ideal 
Eermi gas with equivalent pressures in all directions, though 
the existence of Landau levels. The equation of Pg in a su¬ 
perhigh magnetic field is consequently given by 

1 2 ^ 2 , 

jinpedpe 


Be = 


3/l3 


/o {pIc^ + 


= 1.412 X 10^^^(a;e)dynes cm ^ (12) 

where Ag = is the electron Compton wavelength, Xe = 
ppMI ~ E-p(e) _ gg yy ^ TJL.^ r^Qr )^), and (j>{xe) is the 

polynomial (l>{xe) = g^[a;g(l + Xg)i(^ ~ 1) + + 

(1 + ^^)^]]. 

When p > lO' g cm“3^ Xe 1, and (j){xe) —>■ 

Thus, Eq.(l 1) can be rewritten as 

:) dyneem”^. (13) 


Pg ~ 6.266 X 10^°(-^ ^ 


'po Bee 0.0535' 

It is obvious that Pg increases sharply with increasing P 
when the values of Yg and p are given. 
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3 Magnetic effects on EoSs 

3.1 Magnetic effects on the EoS of BPS model 

By introducing the lattice energy, Baym, Pethick & 
Sutherland (1971) (hereafter “BPS model”) improved on 
Salpeter’s treatment (Salpeter 1961), and described the nu¬ 
clear composition and EoS for catalyzed matter in complete 
thermodynamic equilibrium below pd- BPS model is one of 
most successful models describing matter of the outer crust. 
According to BPS model, the matter energy density is given 
by 


(Shapiro & Teukolskysee 1983); employing /3- equilib¬ 
rium and charge neutrality gives Up = Ue = 9.6 x 
cm“^; /3-equilibrium implies energy conserva¬ 
tion and momentum conservation (pf(p) = PF(e)), we get 
£;F(e) = Pn = E'-pin) = p'^{n)l2mn = MeV, 

and pp = Ep{p) = p|(p)/2TOp = MeV; the 

isotropic matter pressure P is given by 

P = Pe + Pp + Pn 

, , rUpC^ rrinC^ 

= -^nxe) + -^^nxp) -f -^^nxn) , (17) 

Ag Ap A^ 


e = nAr(VFiv(A, Z) + eL{Z, Ue) + ee(ne) , (14) 

where un is the number density of nuclei, Wn{A,Z) is 
the mass- energy per nucleus (including the rest mass of Z 
electrons and A nucleons); Eg is the free electron energy 
including the rest mass of electrons in a unit volume; El is 
the bcc Coulomb lattice energy per nucleus, 

El = -lAUZ^/^e^e^ni^^ , (15) 

where the relations of njv = tibIA and rig = Zbn are 
used. The matter pressure p of the system is given by 

P = Pg + PL (16) 

For a magnetic field B* ':$> 1, Pg in Eq.(15) is given 
by Eq.(12). Based on the above equations, we plot one 
schematic diagrams of QED effects on the EOS of BPS 
model, as shown in Fig. 1. 



Fig.l P vs. p of BPS model below neutron drip. 


3.2 The QED effects on the EOS of ideal npe gas 

We consider a homogenous ideal npe gas under /3- equi¬ 
librium, and adopt ST83 approximation (Shapiro & Teukol¬ 
skysee 1983) corresponding to the weak-field limit as the 
main method to treat EoS of this system in the density range 
of 0.5 ^ 2.5po where electrons are relativistic, neutrons and 
protons are non-relativistic. When neutron pressure dom¬ 
inates, p Ri ninUn, then = 1.7 x 10^®(-^) cm“^ 


where Xp 


pp(p) ^2mp/ip 

rripC rripC 



, the expression 


of (f>{xp) is completely similar to that of (j){xn)- 


Based on the above results, we gain the following useful 
formulae: 


Pp = 1.169 X dynes cm-2 , 

Po 

Pg = 1.825 X 103i(-^)t dynes cm-2 , 

PO 

Pn = 6.807 X 1033(-^)t dynes , 

Po 

Ye = ——— ~ — = 0.005647(-^) . (18) 

Up -\- Un Un Po 


Our methods to treat EOS of an ideal npe gas (sys¬ 
tem) under /3-equilibrium in superhigh MFs are intro¬ 
duced as follows: Combining Eq.(9) with momentum con¬ 
servation gives the chemical potential pp = Ep{p) = 
1.005(-^ — „ MeV, and the non-dimensional vari- 

able Xp = \ ~ 4.626 x 10~^(-^-^ „ )^; Then 

P y rripc-^ 'Bor po 0.0535-' ’ 

we get 





(2 X (43.44( 


B p Ye 
Per po 0.0535 


)l/4 


-1.29-f 1.005( 


B p Ye 
Per Po 0.0535 


)l/2))l/2_ 


(19) 


The /3-equilibrium condition gives the expression for the 
isotropic matter pressure P, 


P = 


■(P[Xe) + —^nXp) + -Ti^(p{Xn) 


A3 


A3 

p B Ye 


A3 


= 6.266 X 10^°(-^ ——+ 2.324 x 10^® 
Po -Ocr U.UooO 

x{ — — -1.624 X 1033 


Po Per 0.0535 


1 


1575^ 


{xl - —xl + —xl) dyne cm 


14 


24 


( 20 ) 


where Xn is determined by Eq.(19). The above equation al¬ 
ways approximately hold in an ideal npe gas when P* ^ 1 
and p ^ 0.5po ~ ‘^Po- Based on the above equations, we 
plot two schematic diagrams of QED effects on EOS of this 
npe gas, as shown in Figs.2-3. Both Pp and P„ increase 
obviously with p and P for an ideal npe gas. 
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Fig. 2 Pp VS. p for an ideal npe gas in a superhigh MR 


The relation of Pp and density p in ideal n-p-e gas 

... - 

= 100 ^ , - ' " ' 



^ " (2). B* « 1 

a.p~(1.4xlO‘'‘-5.6xlO“)g/cm’ ; 


b. Y, = 0.005647x(p/po) 


Fig.5-1 
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oF 
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The relation of and density p in ideal n-p-e gas 


--( 1 ). 5*=100 

^ . - ' " ' a. p ~(1.4xl0‘^-5.6 xlO^^) g/cm^ 
b. Y, = 0.005647x(p/po) 

Fig.5-2 


' • (2). S* 1 
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Fig. 3 Pn VS. p for an ideal npe gas in a superhigh MF. 


3.3 The QED effects on the total matter pressure and 
total energy density 


As discussed above, the pressures of fermions increase with 
B, the total matter pressure increases with B. Due to a pos¬ 
itive co-relation between the total energy density e and the 
total matter pressure, e also increases with B. 

The stable configurations of a NS can be obtained from 
the well-known hydrostatic equilibrium equations of Tol- 
man, Oppenheimer and Volkov (TOV) for the pressure P{r) 
and the enclosed mass m(r). 


dP{r) 

dr 

dm{r) 

dr 


G{m{r) -\- Aiir^P{r)/c^){p -f P{r)/c^) 
r{r — 2Gm(r)lc^) 

Anpr^ , ( 21 ) 


where G is the gravitational constant. For a chosen central 
value of p, the numerical integration of Eq.(21) provides the 
mass-radius relation. In Eq.(21), the pressure P{r) is the 
gravitational collapse pressure, and always be balanced by 
the total dynamics pressure, P; the central density p is pro¬ 
portional to the matter energy density e; the enclosed mass, 
m{r), increases with the central density p when r is given. 

As we know, the magnetic effects can give rise to an 
anisotropy of the total pressure of the system to become 
anisotropic (Bocquet et al. 1995, Paulucci et al. 2011). The 


total energy momentum tensor due to both matter and mag¬ 
netic field is to be given by 

+ T^’', (22) 

where. 


= ernuf^u’'-P^{g^’'-uV), 


and 


B 


1 


rr =- 


B^^B'' 


(23) 


(24) 


47r 2 47r 

The first term in Eq.(24) is equivalent to magnetic pressure, 
while the second term causes the magnetic tension. Due to 
an excess negative pressure or tension along the direction to 
the magnetic field, the component of along the held, 
is negative. Thus, the total pressure in the parallel di¬ 
rection to MF can be written as 


Pn = Pm — 


Stt ’ 


(25) 


(26) 


and that perpendicular to MF, P±, is written as 

n2 

P± = Pm + - A4B , 

Stt 

where A4 is the magnetization of the system, and A4B 
is the magnetization pressure (Perez et al. 2008, Ferrer et. 
al. 2010). In this work magnetars universally have typical 
dipole MFs ^ (10^^ — 10^®) G and inner held strengths not 
more than 10^^ G, under which the system magnetic mo¬ 
ment satishes Ai < B, a condition that can be justihed for 
any medium that is not ferromagnetic, the effect of AMMs 
of nucleons on the EOS are insignihcant and thus ignored 
(Ferrer et. al. 2015). It’s obvious that the total pressure of 
the system becomes anisotropic, that is P± > P||, which 
could lead to the Earth-like oblatening effect. 

According to our calculations, when B* = 100, Pm ^ 
1033 — 1034 dynes cm“^ and ^ ^ 10^® — 10^°. Hence, 
in this presentation, we consider that the component of the 
total energy momentum tensor along the symmetry axis be¬ 
comes positive, > 0, since the total matter pressure 
increases more rapidly than the magnetic pressure. We pro¬ 
pose that the component of the total energy momentum ten¬ 
sor along the symmetry axis becomes positive, since Pm 
always grows more rapidly than the magnetic pressure. The 
magnetic tension along the direction to the magnetic field 
will be responsible for deforming a magnetar along ME, and 
turns the star into a kind of oblate spheroid. Be note that 
such a deformation in shape might even render a more com¬ 
pact magnetar endowed with canonical strong surface helds 
B ^ G. Also, such a deformed magnetar could 

have a more massive mass because of the positive contribu¬ 
tion of the magnetic field energy to EOSs of a magnetar. 


4 To modify Pe in superhigh MFs 

According to atomic physics physics, the higher the orbit 
quantum number I is, the larger the probability of an elec¬ 
tron’s transition (this transition is referred to the transition 
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from a higher energy level into a lower energy level) is. 
Analogous to atomic energy level, in a superhigh MF, the 
easier an electron’s transition from a higher Landau level 
into a lower Landau level. Thus, the higher the electron 
Landau level number n, the lower the stability of the Lan¬ 
dau level. Owing to the uncertainties of microscopic states, 
we introduce a new quantity, 5 „, the stability coefficient of 
electron Landau level in a superhigh MF, and assume that 
Qn decreases with n as an exponential form, 

5n = <?on“, (27) 

where po is the stability coefficient of the ground-state Lan¬ 
dau level of electrons, a is the Landau level stability index, 
and is restricted to be a < 0. From Eq.(27), it is obvious 
that g{n) is a function of n and a, and the higher n is, the 
smaller is (except for gi = go)- 

According to the Pauli exclusion principle, electron en¬ 
ergy state number in a unit volume, Npi^a, should be equal 
to electron number in a unit volume, rig. Considering the 
electron Landau level stability coefficient g^, and summing 
over electron energy states in a 6-dimension phase space, 
we can express Npha as follows: 

^pha — — ^AP^e 



X 6{— - [{2n + 1 + a)B*]^)p±dp±, (28) 

J rUeC 

where Na is the Avogadro constant. When rim > 6, the 
summation formula can be approximately replaced by the 
following integral equation 


E 

n—O 


n“+= ~ 


n 


“+^dn = 




2o; -f 3 


(29) 


Thus, Eq.(28) can be rewritten as 

Npha = NApYa = 


r —— , TTlaC . q 

7r/^(—— 


fEF{e)^2 1 / Pz P 


rrieC 


d(^). (30) 

rrieC 


Jo rrieC" 

After a complicate deduction process, we get an non- 
dimensional momentum 


PF{e) 

rjif-c 


= C 


Ya 


2(c + 2) 


(31) 


0.05 po 

where C is a constant, which is determined by 

0.05po-^A(2a + 3) i h _a 

C = (- , , -i) 2{c + 2) { -) 2(c + 2) 

22(i-“)7rpo/(a) ^rrieC 


~ (337.12)5(^(- 


2ct -f 3 


-) 2(c + 2) _ 


’22(i-“)go/(a)' ' 
with /q^(1 — ^ _ p^c/Ep{e). If a and go are 

determined, the expressions of EF{e) and Pe in superhigh 
MEs will be modified accordingly. To exactly determine the 
values of a and go is an interesting and important task, but is 
beyond of this paper. Since this is an original work in which 
some uncertainties could exist, to further modify and perfect 
our model should be considered in our future studies, espe¬ 
cially to further investigate QED effects on the EoSs using 
an improved expression of Pg in a superhigh ME. 


5 Conclusions 

In this presentation, we derived a general expression for 
electron pressure, which holds in a superhigh MF, consid¬ 
ered QED effects on EoSs of neutron star matter, and dis¬ 
cussed an anisotropy of the total pressure in superhigh MEs. 
Compared with a common pulsar, a magnetar could be a 
more compact oblate spheroid-like deformed NS, due to the 
anisotropic total pressure; an increase in the maximum mass 
of a magnetar is expected because of the positive contribu¬ 
tion of the magnetic held energy to EoS. 
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